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Abstract Interpolation by various types of splines is the standard procedure in many 
applications. In this paper we shall discuss harmonic spline "interpolation" (on the 
lines of a grid) as an alternative to polynomial spline interpolation (at vertices of a 
grid). We will discuss some advantages and drawbacks of this approach and present 
the asymptotics of the Lp-error for adaptive approximation by harmonic splines. 
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1 Harmonic splines as an alternative to polynomial splines 

In Is) authors investigated the question of adaptive approximation by polynomial 
splines on box partitions in W'. The obtained a general description as well as sharp 
constant in several popular cases. However, polynomial splines might not be the best 
option when it comes to interpolating a function over box partitions. There is an 
alternative - harmonic splines - which can be viewed as a direct multivariate gener- 
aUzation of linear splines in one dimension. This paper addresses approximation of 
bivariate functions and, therefore, from now on we shall consider harmonic splines 
of two variables. 

The following several observations help to see the analogy between univariate 
linear splines and bivariate harmonic splines on rectangular partitions more clear: 
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1. Linear univariate splines interpolate the given function on the whole boundary 
of the mesh which is the points of the partition of the domain. While bilinear 
splines defined on rectangular partitions interpolate the given function only at 
the vertices of the mesh elements, not on the whole boundary, harmonic splines 
are constructed to coincide with the original function on the boundary of the 
meshes - on the lines (as oppose to just at the points). Litvin in [10] introduces 
different terminology to distinguish between these types of interpolation and he 
says that harmonic splines "inter/Zneafei" function instead of just interpoZaring. 
Using harmonic splines has advantages and drawbacks and the choice depends 
on a particular problem. The strong advantage is that it uses more information 
about the function. However, it also imposes restrictions on the mesh. In the case 
of square domain (considered in this paper) the meshes will consist of mostly 
squares. In the case of a rectangular domain, the elements of the mesh will have 
the same shape (same width to length ratio) as the domain, and only size will 
adapt to the local behavior of the function. 

2. Both univariate linear spline Pi{x) and bivariate harmonic spline U{x,y) satisfy 
corresponding operator equations. Indeed, if we denote by Zii = ^ and A2 = 
■^ + ^, then clearly 



3. A univariate linear interpolating spline and bivariate harmonic spline attain max- 
imum and minimum values on the boundary of the element. 

4. If the values of the univariate linear spline at internal nodes are determined from 
the continuity condition on the first derivative, spline degenerates into a segment 
of straight line. If the values of harmonic spline are determined from the condition 
of continuity of the normal derivative to the interior boundary of the mesh, then 
the spline degenerates into a harmonic function on the union of domains (which 
is, by the way, far less restrictive). 

5. Both have very similar local estimates of the error of approximation JSj. 

For more discussion of properties of harmonic and polyharmonic splines (or 
polysplines) see the monographs ll9l [T0l and references therein. 

To approximate the bivariate function (so that the approximant also coincides 
with the given function on some fines) polynomial splines, in particular, blending 
splines, can also be used. However, tfiere exists a broad class of functions tfiat can 
be well approximated by fiarmonic splines, but are not at all approximated by such 
polynomial splines under the same requirements concerning the order of smoothness. 
Let us present an example of such a function. 

Example ( 13). Let us consider the function 



AiPi{x) = and A2U{x,y) = 0. 




{smh{ak{y - b)) + sinh{akb) 



smh{aky)) sin(a<.jc) 



0<x<a,0<y<b, {2k+l)K/a. 
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For this function 

lun = oo, 

x->0,y-!-0 dxay 

and, hence, it cannot be approximated by blending polynomial splines. On the other 
hand, Af(x,y) = — 1 which implies that / can be well approximated by harmonic 
splines. 

In univariate case the idea of using splines constructed based on solution of some 
differential operator equations, called L-splines, has been more or less well-developed 
(see, for example, ifTSll and references therein). There has been several attempts to 
generalize this concept to multivariate case. For example, Litivn in his monograph 
IfTOl introduces concepts of "interlineation" and "interflatation". These concepts are 
natural and direct generalizations of interpolation to the case when one needs to re- 
construct the given function based on knowing values of function (and possibly its 
normal derivatives up to some order) on one or several lines ("interlineation") or m- 
dimensional (1 < m< d) linear manifolds ("interflatation") in E''. 

As already mentioned earlier, in IfTOll for the given on function Litvin sug- 
gests and compares two methods of polynomial and non-polynomial "interlineation" 
on several lines. The first (polynomial) is a generalization of Hermite polynomials 
which preserves the smoothness order of the original function. The second (not poly- 
nomial) operator is a generalization of D' Alambert operator, which is the solution of 
homogeneous wave equation, to the case when data is given as values of the normal 
derivatives of aU orders up to and including N {N > I). 

There exist many applications in which data is given not at a discrete set of points 
but on some linear manifolds (lines, planes, etc.). One of many examples is cartogra- 
phy. The measurement and analysis of bathymetric measurements is one of the core 
areas of modern hydrography, and a fundamental component in ensuring the safety 
of water transportation, efficiency of offshore resource development, precision and 
effectiveness in recovering underwater objects as well as search and rescue missions. 
The data used to make bathymetric maps today typically comes from an echosounder 
(sonar) mounted beneath or over the side of a boat, "pinging" a beam of sound down- 
ward at the seafloor or from remote sensing LIDAR or LADAR systems. In all the 
cases, the data obtained are univariate - along the lines, curves or, more generally, a 
grid composed of them - depending on the course of the surveyor vessel (which can 
be adapted when needed). The goal is to recover the function describing the relief of 
the seafloor based on the "traces" of this function on the lines which described the 
course of the vessel with sonar 

Similar technology is used for mapping surfaces of cosmic objects. Only in- 
stead of sonar, which uses sound waves, radars are used. Radars are based on high- 
frequency electromagnetic radiation which can penetrate the surroundings of objects, 
for instance, the clouds covering Venus. 

Another example of applications to surface reconstruction based on data given 
on lines, curves, or hyperplanes would be manufacturing of car, aircraft etc. bodies. 
Currently, most popular tool used in these models is polynomial splines. Being easy 
to manipulate, they have many drawbacks: lack of global smoothness (if working with 
convenient low degree polynomial pieces) and lack of flexibility in choosing the mesh 
(for instance, it might be necessary to base a part of the body on hexagon or pentagon. 
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not only triangle or rectangle). Harmonic and their generalization - poly harmonic - 
splines have great advantage in both directions: they preserve higher smoothness and 
they can be constructed on any (more or less reasonable) domain. 



2 Notation, statement of the main problem, and the main result 

Let in the space K.^ of points x = (xi ,X2) the unit square [0, 1] x [0, 1], be given with 
interior denoted by I^, and the closure denoted by F. By Dn we shall denote a par- 
tition of P whose number of elements has order N as N ^ °°, and is so that the 
majority of elements are squares, however some small number {o{N) as N ^ °°) of 
rectangles is also allowed. The interior of an arbitrary element of the partition On we 
shall denote by f2^, the closure by £2 , and the boundary by dQ^ . 

We shall need the following standard notation: C{£2) denotes the space of contin- 
uous in some region £2 functions, C^{Q ) denotes the space of twice differentiable in 
Q functions, and Lp{Q), I < p < °° denotes the space of measurable and integrable 
in power p functions f : Q ^M. with the norm 



In order to introduce the set of functions that we will use as an approximation tool, 
let us denote by H{Q) the set of harmonic functions in bounded domain £2 C M?, i.e. 

H{Q) = {m(x) e C^{Q) : Au{x) = 0} , 

where A = Jjr + is the Laplace operator 

Now we shall consider the set 5'(nAf) of continuous on P functions such that 
when restricted to any 12^ from the partition On are harmonic functions, i.e. 

SiON) = {ge C(F) : g\^^ = u.l^N, ui G H{Q'') Vr2^ G D^v} • 

From now on we shall refer to the functions from set 5(nAr) as harmonic splines. 

Let now a function / € C^{¥) is given. We will approximate it by a harmonic 
spline s{f, On) in such a way that their values coincide on the boundary of elements 
^2^^. In other words we shall require the spline s{f,ON) to "interlinate" the given 
function on the boundary of all partition elements. Hence, the function s{f, On) in 
each shall satisfy the equation 

As{f,aN)^o 

and the boundary condition 

In other words, function s{f, On) in each domain On must be a solution of Dirichlet 
problem for Laplace equation, which is unique (see, for example, ifTSl ) and therefore 
harmonic spline s{f,ON) for each / e C^(P) is well defined on P. 
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The main goal of this paper is to study the minimal Lp-error (1 < p < °°) of 
approximation of a given function / by harmonic splines s{f,\I\N) for all possible 
partitions Di\i which we shall denote by 

RNif,Lp)=mf\\f-s{f,nN)h,^^py 

The value of RN{f,Lp) we shall call the optimal error. The partition on which 
the optimal error is achieved will be called optimal partition, i.e. 

/?^(/,L,) = !i/-.(/,a;v)iL^(i2). 

For an arbitrary function / it is impossible (except for some trivial cases) to find the 
optimal partition 'OiN and explicitly compute the value of the optimal error. Therefore, 
we are interested in the following two natural questions: first of all, how the optimal 
error behaves as number of partition elements increases (A^ — )• oo), and secondly, how 
to construct a sequence of partitions {D^l^^j, that is asymptotically optimal, i.e. 

||/-.(/,n^)||^(i2) 

lim — i ~ 1. 

RNif,Lp) 

In order to state the main result of this paper, we need the following notation. By 
Gq (x; v) with x, \ E £2 we shall denote the Green's function of the Dirichlet problem 
for the domain Q (for the detailed definition, see, for example, ifTSl ). In addition, 
denote by 

/(x) := J Gj2(x;v)c/v, x e I^. (2.1) 

The main result of this paper is the following theorem. 

Theorem 2.1 For an arbitrary function f e C^(P) there exists a sequence {Dj^UJ^j 
and corresponding sequence of splines such that 




Remark 1. Note that the majority of the partition elements from D)^ are squares. 
This is rather restrictive condition on the mesh. However, it is somewhat natural to 
expect since we are using more information about the function than in the case of 
interpolation by polynomial splines. 

Remark 2. In the case when the domain of the function is rectangle, the mesh 
elements must have the same shape as the domain, i.e. the same width to length ratio. 

Comparing this result with analogous results in the case of approximating the 
given function from C^(P) by bilinear polynomial splines (|[3]|5l), we see that the or- 
der of the error of approximation is the same. Therefore, depending on the particular 
problem one could choose either polynomial or harmonic splines, the latter providing 
interpolation on a larger set, but allowing less flexible partitions. 
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3 Idea of the proof and auxiliary results 

The main idea of the construction of each partition in the sequence is to use first "in- 
termediate approximation" of the given function / which needs to be approximated 
by splines. As an intermediate approximation we shall use the piecewise function 
"glued" from second degree Taylor polynomials of /. Then the intermediate approx- 
imation is used to build the asymptotically optimal sequence of partitions of the do- 
main and, consequently, the sequence of corresponding harmonic splines. 

The process of constructing this sequence of partitions consists of two steps: 

1 . The domain P is first divided into some number (small comparing to A^) of equal 
sizesubdomains, and on each instead of the original function / we consider its 
second degree Taylor polynomial constructed at, say, center of the corresponding 
subdomain. The number of these subdomains is chosen so that the absolute value 
of the difference between the function and its Taylor polynomial is small enough 
on each corresponding subdomain. 

2. Next, each subdomain is refined further, depending on the Taylor polynomial 
(and therefore the original function) on this subdomain. The final partition D'^ 
will consist of the squares and possibly some rectangles (their number is small 
comparing to A'). The total number of partition elements is determined from the 
following two conditions: first of all, the total number of all the elements in D)^ 
is approximately equal to A^, and secondly, the global error of approximation of 
piecewise quadratic function (consisting of Taylor polynomials built on the pre- 
vious step), by corresponding harmonic splines is minimal. 

The problem of finding this error of approximation on "smallest' subdomain ele- 
ments is reduced to approximating the function Q{x) = Ax\ + Bx^, where A and 
B are constants, since the rest of the terms in second degree Taylor polynomial 
constitute a harmonic function. 

The constructed in such a way partition for each fixed is used in the proof 
of the main Theorem l2.1l 

The idea of this construction based on intermediate approximation by piecewise 
quadratic functions has been already used in papers HI, ||2l, Q, Q by authors 
and co-authors to build asymptotically optimal sequences of polynomial splines on 
triangulations, rectangular partitions, and their generalizations in various settings. 
More on the history of adaptive approximation and asymptotically optimal sequences 
of splines in the case of polynomial splines can be found in IJl, [S). 

Let us turn now to the auxiliary results that we shall need to prove the main 
theorem of this paper. 

For a function / e C2(P) we shall define a modulus of continuity as follows 

a)(/,5)=sup{|/(x)-/(x')|: ki < 5, |x2 -x'jj < 5, x,x'gF}. 
Then we consider 



to ( 5 ) = max { 0) (/xv , 5 ) , 0) {fyy , 5 ) , o) (/vj , 5 ) } 



(3.1) 
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Lemma 3.1 Let function f G C^(P) be given and let -P2(x) = P2(/;x;xo) be its sec- 
ond order Taylor polynomial taken at the point xq, that is the center of some square 
Df, C P with side length h. Then 

|/(x)-^2(x)|<^«f^y xSA,. (3.2) 



2 V2^ 

For the proof of this lemma see, for example, IT]. 

To state the next lemma we shall need the following notation. 

Let Gq (x; v) be a Green's function of the interior Dirichlet problem for bounded 
region i2 C M^. By 12 + d we denote the region obtained by the shift of the region 
^2 by a vector d E M", and by aQ we denote the region obtained by scaling £2 by 
coefficient a. 

Lemma 3.2 The following properties of the Green's function hold true: 

Gaf2(x;v) = Gn ( — ; — ) , x e a£2,\ e aQ, (3.3) 
\a aJ 



Gi2+rf(x;v) = Gi2(x-d;v-d), xe f2+d,ve f2+d. (3.4) 

Proof. In the proof of this lemma we shall use the well-known (see, for in- 
stance, lfT3]| ) formula for the Green's function of the interior Dirichlet problem for 
the bounded region Q . For xE Q, \ E Q the Green's function can be written as 

^^("'") = ^^"h^+^^"'")' 

where function g{x;\) is harmonic in both arguments on Q, continuous in v on 12, 
and is chosen so that the Green's function has zero value on the boundary. 
Forx e aQ, v e aQ 

xv\ 1. 1 /xv\ 1. a /XV 



a al 2k - — - \a a) 2k |x — v| \a a 

= -i- Mna + ln-^^ — - ) +g {—■,—\ = ^\n-—^ — -+|(x;v), 
27r V |x-v|/ \a a) 2k |x-v| 

where function g(x; v) = ^ (§; ^) + In a is clearly harmonic in both arguments in 
aQ and continuous in v on aQ . Besides that, since the boundary of the region Q is 
mapped to the boundary of the region aQ, and the Green's function has zero value 
on dQ, then the expression above will be zero on d{aQ). Hence, ( 13.31 ) is proved. 
Let us consider next the shift of the region Q . For xei2+d, vei2+d 

G^(x-d;v-d) = i^ln^— ^^1^— ^+,(x-d;v-d) = 

= :^lni x;v). 

2k X — V 



8 



Yuliya Babenko, Tatyana Leskevich 



Since all the necessary conditions for the function ^(x; v) = §(x — d; v — d) are satis- 
fied, then the obtained expression is the Green's function for the region £2 +d, and 
hence ( 13.41 ) is proved. □ 

Next we shall consider a square region f2 C M" with sides parallel to the coordi- 
nate axis, and on which the following functions is defined 

Qix)=Ax] + Bxl (3.5) 

where A and B are constants. 

We shall find the Lp-error of approximation of the function Q{x) by a harmonic 
function u{Q;x) on Q such that the values of u{Q;x) coincide with the function 
Q{x) on the boundary of Q. In other words, for the function u{Q;x) the following 
conditions must be satisfied 

Au{Q;x)^0,xen (3.6) 

and 

uiQ;x)\dn^Qi^)\dn- (3.7) 

Lemma 3.3 For the quadratic function Qix) = Ax\ +8x2 such that AQ{x) = 2(A + 
B) 7^ 0, we have 

llG-«(G)llMi2)=2|A + BMi2|'+i||/||^^^(,2), (3.8) 

where \Q \ denotes the area of Q, and I is defined in ( 12.71 ). 

Proof. We shall use the following representation of an arbitrary function «(x) € 
C^(i2) nC(i2) on the region i2 C with smooth enough boundary: 

dGn (x;v) 
da 

da n 



m(x) = - / — ' u{\)d\~ I Gn(x-\)Au{\)d\,x€Q, (3.9) 



where '^'^"jf''"'* is its partial derivative in the outer normal direction to dQ (see, for 
example, ||T31 ). 

Then for the function Q{x), taking into account i3.9i and the fact that AQ{x) = 
2{A + B), we have 

Q{x) = - j ^^I^^Q{y)dy-2{A+B) J Gn{x;\)dy, x e Q. (3.10) 

dQ Q 

In addition, with the help of (|3.9l l, the solution of the problem (I3.6l l - ( I3.7l i can be 
written as 

dOo (x;v 



u{Q-x)^- ^^e(v)fl!v, xei2. (3.11) 

J dn 

dn 

From (13.1 Ok ( 13.1 lb . and the fact that the Green's function is nonnegative, it fol- 
lows that 

|e(x) - u{Q;x)\ = 2\A+B\ j Gn{x-\)dy, xeQ, 
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and, hence, 



\Q-u{Q)\\L,(n) 




GQ(x;\)d\ dx 



(3.12) 



Our next goal is to transform the integral in the right-hand side of ( I3.12l i so that it ex- 
plicitly depends on the area of i2 . In order to do so we need the result of Lemma l372l 
Let f be the square unit region such that when centered at one of its vertices 
and stretched by a it becomes Q. Then the volumes of these regions are related by 



1X21 



a jl"!, from where we have a — v I'^l- 
Therefore, we have the following equality 




Gn(x\\)d\ dx 
\n \n I 

Changing variables we arrive at 




p 



GQ{x;\)d\ dx 



G,J^;^)dy 



dx 



\n \n 




Gy2(x;v)c/v dx 



(3.13) 



.12 \l2 



It is left to notice that due to property ( 13.41 ) of the Green's function, the last equal- 
ity will not change after the shift of the region by an arbitrary vector from M?. 
Therefore, as we may take an arbitrary rectangular region of unit area with sides 
parallel to the coordinate axis, in particular, we may take I^. 

Taking into account ( 13.131 ). notation ( 12. Il l, and having recalled that for the given 
region Q the value a = \/\^ |, equality ( I3.12l i becomes 



lG-«(e)ilMi^)=2|A + BMi2|'+^||/|L^^(i: 



□ 



Now let us turn to the proof of Theorem 12.1 



4 Proof of the main result 

We shall begin the proof by constructing (for a given function / e C^(P) and fixed 
number A^) the partition of P such that the sequence {D^l^^j is asymptotically 
optimal. The idea of this construction is given in Section[3] 

First of all, we determine the number of the elements of the intermediate partition. 
To this end, for an arbitrary fixed £ G (0, 1) we set 
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where (o{5) was defined in ( 13.1b . It is clear that °° asN 



We shall take as the number of elements of intermediate partition, and we 



shall subdivide the square P into equal squares Djy , / = 1 , . . . , with the side length 
. Let us show that the number of elements in intermediate partition is indeed small 
comparing to the total number of elements A^, i.e. 

mjj = o{N), N^oo. (4.1) 

In order to do so let us consider : 

»'iv 

1 / 1 \ 2(mN-l) 



I, 2(m^-l)2 \2imN-l)J fn2(^ 

N I2(mjv-1) 



>e ^-^ (4.2) 



mlco 



Taking into account ( '-^^^lA ) _^ i and (O ( ) ^ as A^ ^ oo we obtain -^^! — > o° 

as A^ — > oo, and hence (14.1b holds true. 

In order to obtain the intermediate estimates we define the functions /Af(x) and 
Qn{^), using the notation (/; x; x/ ) for Taylor polynomial of second degree for /(x) 
at the center x/ of each square D^, as follows: 



1. ForxeD^set 



and 



/w(x)=P2(/;x;xi) 



2. When 1 < / < m^, for x G Z)f \ u'-ZldOf we set 



/^(x)=P2(/;x;x,) 



and 



Next we shall find the number of elements to additionally subdivide each Djy, I = 
l,...,m^. For that we shall cover the square Djy by a "mesh" consisting of squares 
of the same fixed area, fixing, for convenience, one of the mesh vertices with one of 
the vertices of the square Djy. The intersection of Djy with this mesh will provide the 
necessary subdivision of the square. It will consist of the squares of the original mesh 
as well as, possibly, some rectangles along the boundary of Djy. 

We shall consider two cases: when the partition consists only of squares and when 
the partition contains some rectangles along the boundary. 

First, let us assume that partition of Djy consists of squares only. We shall find the 
area of the squares next. The number of the squares in the partition of Djy we shall 
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denote by hj. The elements itself we shall denote by Q'^', i ~ l,...,fif, I ~ 1, 



^2 



Then |J |J Q^' gives the partition of (for fixed A^), and hence =N. 



I=\i=\ 



1=1 



To determine nj for each D{y we shall minimize the value of the error of approx- 
imation of /Af(x) by the corresponding harmonic spline s{fN, D]^) on P. The needed 
error can be written in terms of the errors on each element of the partition as 



1=1 i=l 



(4.3) 



Therefore, we need to find the error of approximation of /Ar(x) on each element X2^'. 
Since on 12^' the difference /Af(x) — Qat (x) is a harmonic function, using Lemma [331 
we have 



wiw: 



where we used the notation 



Taking into account ( 14.41 ). the equality ( 14.3b becomes 



(4.4) 



(4.5) 



II/. ^ s{fN, D^)!!^^,,.) ^ LKNi-') (^,,;)2(;,.l) \K„i,^f> 

from where we finally arrive at 

...2 

1 



/;v-.K/.,n^)iij;^^(i2)-EM;j,(x,)^^^^2(p+i)(,,,)2;,ii'iiMP)- 



(4.6) 



Using the method of Lagrange multiphers, we shall minimize the right-hand side of 
(I4.6l l under constraint 



(4.7) 



1=1 



For that let us consider the function 



^(ni,..,n,,,2,X) = Y^M'l^ix,)j—^ 
where A is the Lagrange multipher. 
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To find the critical points of the function 0{n\, ..,h^^2 ,A) we shall consider the 
system 



„2 



/=i 

Solving it we obtain hj: 



(4.9) 



^'"^^'\/ = l,...,,4- (4.10) 



im;;'(x,-) 

1=1 ■' 

Since for the solutions of the system ( 14.91 ) the value d^^" > 0, then ( 14.1 Oi l indeed 
provide the minimal value of the total error of approximation of the function /Af(x) 
by harmonic spline i(/Ar, D]^) on the square F. 

As the area of every square of the mesh Q. C Z)jy for the fixed we shall take 

[niNnir 

Later we shall need the following estimate for the values «/ 



V^min{M;|-"(x)} ^ 
ni> — =C^>0, (4.12) 

II A rll 2(p+l) 'Wa? 

where C is independent of constant. 

For the final partition of the domain to consist of only squares, the values defined 
by (14.10b for all / = 1, ...,ot^ must be squares of natural numbers. In this case we will 
have optimal partition. Otherwise, the final partition will contain both squares and 
rectangles, and will be optimal only asymptotically. 

From now on we shall assume that the partition of Djy consists of both squares of 
the mesh and some rectangles along the boundary. By we shall denote the set of 
all squares, and by - the set of rectangles from the partition of Djy. In addition, let 

~ "?(^) denote the number of squares from Then the number of rectangles 
from Kj^ is clearly equal to 2«/ + 1 . 

It is also clear that for all / = 1, ...,m^ we have 

{n,-\f <n1 <n]. (4.13) 

The number of elements of the partition of D'^ (in the case when rectangles are 
present) is («/ + 1)^. Let us show that the total number of the partition elements 
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+ and are values of the same order asN To that end, let us estimate 

1=1 

1 2 

the values nj and «/■ Using inequality (14.131 1 and condition ( 14.71 ) we have 

1=1 i=\ 



ml 


ml 






1=1 


1=1 






ml 


ml 






1=1 


1=1 



2 2 2 2 

nf, mj, mf, mj, 

5^ „2 



(4.14) 



(4.15) 



1=1 1=1 1=1 1=1 
Using Holder inequality, we obtain 

ml ml I ml \ ^ 

^ «/ < ^ «/ < niN L "M = "^nVn. (4. 16) 

/=i i={ \i=i I 

Taking into account estimates ( 14.141 ) - ( 14.16b . we have the following double inequal- 

ml 

N - IiiinVn < XI («/ + 1 )" < ^ + ^ihnVn + ml. 

1=1 

Dividing both sides by A^, taking the limit as A' ^ oo and using ( 14.11 ). we obtain 

ml 

I («/ + !)' 
lim^^^ = 1. 

Therefore, the number of elements of the constructed partition has order N as N ^ °°. 

Next we shall verify that the constructed in such a way sequence of partitions will 
be asymptotically optimal. 

Using the triangle inequality we have 

\\f-s{f,DN)h^,^j2^<\\f-fNh,M^) + \Hf,nM)-s{fy,aN)h,,ii2) 

+ \\fy-s{fy,nN)h^,^py (4.17) 

We shall now estimate each term on the right. 

For the first and second terms we shall first obtain the estimates on an arbitrary 
square Djy of intermediate partition in terms of \\f — /nW^ (q/ y Since — 
then the following is true 



\\f-fN\\'o,<-^\\f^fN\\',y (4.18) 

To estimate the difference D^) — 5(/Af,n^)||^ j^^/ j we shall consider an arbi- 
trary element of the partition Q. C Z)jy. According to the definition of a harmonic 
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spline, the difference of functions s{f, □i^)(x) — s{fN, □^)(x) is a harmonic function 
on Q and its values on dQ coincide with the values of the difference /(x) — /Af(x). 
Therefore, using representation ( 13.9b . we obtain 

D^) (x) -s{fM){x)^-J ^^^^ (/(v) - /^(y) , X e i2 . 

dn 

Then 

K/,n^)(x)-.(/;v,n^)(x)|<||/-MI,.(^) J ^^^^^dy,xen. 

dn 

Taking into account the uniqueness of the solution ( 13. 9t (see, for instance, lfT3l ). we 
have 

J dn 

dn 

and, therefore, for all x e Djy 

\s{f,a*^){x)-s{fy,a*^){x)\ < Wf-fNh^^i^^y 

which implies 

\\s{f,n*^)-s{fy,n*^)\\l , <^\\f-fN\\l (4.19) 

Using now Lemma lTTI and definition of m^, we obtain 

The last inequaUty, together with ( 14.191 ) and ( 14.181 ), imply the needed estimates: 

!l/-/^v|L,,(i2)<^, (4.20) 

||.(/,n^)-^(/^,n^)||^^(i2)<^. (4.21) 

Next let us consider the third term in ( 14.171 1 

ml 

WfN - sify, Dm)\\1^,2) = L WfN - s{fy,aN)\\l^^^,^^ 

"4 / \ 

= L L \\fN-sifN,Dj,)\\l^^^^^+ ^ \\fy-sify,Dj,)r^ \. (4.22) 
1=1 \neRi^ neR'^ " / 

Since for each rectangle Q E R'^^ there exists such an element of the mesh Q that 
Q C £2, then 

ll/^-^(/;v,n;v)li;;^ ^ < |l/^-^(/;v,n;v)||^,,(^). 
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Therefore, using equality ( 14.4b for expression ( I4.22l i we obtain 



Taking into account the estimate ( 14.131 1 for the number of elements we have 

II ^ ff n Ml/' ^ II^I^Mi') I ^("rflli 



ni n. 



We shall estimate nj'' using (14.1 Oi l, and for «/ and nf we shall use ( 14.12b 



1 \ p 



(x,) 

„2 



Since 



\ f m/^ (x,-) ^ / 4/7^ (x)flfx, A? ^ 

,=1 ' 

then for large enough 



\\fy-s{fy,nM)\\l^^,^<^lf^l (1+e) / 4/^^(xyx (^1 



Going back to ( 14.171 ). and taking into account ( 14.201 ) and ( 14.21b . we obtain 

ll/-(/n»)llv-.^|+^((' + ^)/^/*W<'«) ' (> + w4'' 

Hence, using (14.11) . we have 

limsupA^||/-.(/,n^)||i^^(i2) =2£ + ||/|L^^(i2) (1+e) / Af^ix)dx 



Ell 

p 
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Finally, since e e (0, 1 ) is arbitrary, then switching to the limit as e ^ +0, we obtain 

limsupA?||/-5(/,nA,)||, (j2) = PIIl (i2)||4/!Il^(i2). 

In order to obtain the estimate from below we shall use again the function /Af(x) for 
intermediate approximation. We shall use the triangle inequaUty for \\fN — s{fN, 
as follows 



11/ -^(/,n^) 11^^(12) > \\fN-s{fN,n*,,)\\L„ifi)-\\f-fN\\L„ifi) 

-\\s{f,U%)-s{fy,U%)\\^^^^^^y (4.24) 

The estimate from below for the first term can be obtained very similarly to the way 
we obtained the estimate from above for this term. The result of it (for large enough 
A^) will be 

/ \ — 1 

ll/.-(/.,n^)M)>^f(i-.)/l^/(x)l^^xj " (i-^)'. 

Now we may rewrite ( 14.24b . using for second and third terms the inequalities ( 14.201 ) 
and ( 14211 ) 

Multiplying both sides by A^, switching to the limit asN and taking into account 
(lO 

/ \ ^ 

liminfA^||/-.(/,n^)||^^,(i2) > ||/|L,,(i2) f (1-^)/ \m^)\^dxj -2e. 

Since e is arbitrary, we arrive at 

liminfA^||/-,v(/,n^) 11^^(12, > \\Ih„(i2)\\mL_^ii2y 

Together with the estimate from above it completes the proof of Theorem l2.1l □ 
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